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We investigate the multiquantum vortex states in type-II superconductor both in "clean" and 
"dirty" regimes denned by impurity scattering rate. Within quasiclassical approach we calculate 
self-consistently the order parameter distributions and electronic local density of states (LDOS) 
profiles. In the clean case we find the low temperature vortex core anomaly predicted analytically 
in G.E. Volovik, JETP Lett. 58, 455 (1993) and obtain the interference patterns in LDOS distri- 
butions analogous to found previously in the framework of Bogolubov-de Gennes theory. In dirty 
regime the multiquantum vortices feature a peculiar plateau in the zero-energy LDOS profile which 
can be considered as an experimental hallmark of multiquantum vortex formation in mesoscopic 
superconductors. 
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I. INTRODUCTION 

Modern technology development provides a unique 
possibility to study superconducting states at the 
nanoscale. Recently there has been much experimental 
effort focused on the investigation of exotic vortex states 
in mesoscopic superconducting samples of the size of sev- 
eral coherence lengths^. Magnetic field can penetrate 
the sample in the form of a poligonlike vortex molecule 
or individual vortices can merge forming multiquantum 
giant vortex state with a winding number larger than 
unitj«2. The latter possibility is of particular interest 
and the search of giant vortices in mesoscopic super- 
conductors was performed by means of various experi- 
mental techniques including transport measurements^, 
Bitter decoration^, magnetometrjji, and scanning Hall 
probe experiments^. Currently much effort is invested 
to the studies of nanoscale superconducting samples 
with the help of scanning tunneling microscopy (STM) 
techniques'^ which have been achieved only recently 
and allows for the direct probe of the structure of vor- 
tex cores through measurement of the electronic states 
LDOS distribution modified by vortices. 

Such STM measurements have been proven to be an ef- 
fective tool of experimental study of electronic structure 
of vortices in bulk superconductors^— . Indeed for the 
temperatures much lower than the typical energy scale 
in superconductors T <C T c the local differential conduc- 
tance of the contact between STM tip and superconduc- 
tor as a function of voltage V: 
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where (dl/dV)^ is a conductance of the normal metal 
junction and Nq is the electronic density of states at the 
Fermi level. The observation of the zero-bias anomaly of 
tunneling conductance at the center of singly quantized 
vortices^— clearly confirmed the existence of bound 
vortex core states predicted by Caroli, de Gennes and 



Matricon (CdGM}^£. In clean superconductors for each 
individual vortex the energy s(fi) of a subgap electronic 
state varies from — A to +A = A(r = oo) as one 
changes the angular momentum [i defined with respect to 
the vortex axis. At small energies |e| Ao the spectrum 
is a linear function of \i: 



■(//) = U3fl 
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Here u ~ Ao/(fcF£) where £ = HVf/ Aq is coherence 
length, kp is Fermi momentum and Vp is Fermi velocity. 
The wave functions of the subgap states are localized 
inside the vortex core because of the Andreev reflection 
of quasiparticles at the core boundary and determine the 
low energy LDOS singularity at the vortex center. 

In multiquantum vortices the spectrum of electronic 
states bound in the vortex with the winding number M 
contains M anomalous branches degenerate by electronic 
spii 
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where Uj ~ Ao/(/cp£), index j enumerates different spec- 
tral branches (1 < j < M), —kp£ < fj,j < kpS,- Each 
anomalous branch intersects the Fermi level and con- 
tributes to the low-energy LDOS. The spectrum of local- 
ized electronic states in mesoscopic superconductors with 
several vortices have been shown to be very sensitive to 
the mutual vortex position^. It has been suggested that 
testing the properties of electronic spectrum by means of 
the heat conductivity measurement one can directly ob- 
serve the transition to the multiquantum vortex state in 
mesoscopic superconductor—. An alternative route is to 
use STM measurement of local tunnelling conductance 
being proportional to the LDOS provided T <C T c . Thus 
to provide the evidence of multiquantum vortex forma- 
tion revealed by STM experiments one should find dis- 
tinctive features of the order parameter structures and 
LDOS profiles occurring especially in the low tempera- 
ture regime r«T c . 

Previously the low temperature properties of multi- 
quantum vortices have not been investigated much. The 
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results of theoretical studies are known only for the par- 
ticular case of vortices in clean superconductors when the 
electronic mean free path is much larger than the coher- 
ence length. In this regime the contribution of anomalous 
branches produces singularities of the order parameter 
distribution near the vortex core in the limit T « T c . In 
particular the singly quantized vortex features an anoma- 
lous increase of the order parameter slope at the vor- 
tex center which is known as Kramer-Pesch effec t 25 ' 26 . 
The generalization to the multiquantum vortex case was 
suggested in Ref. (f27h where it was analytically predicted 
that doubly quantized vortex should have square root sin- 
gularity of the order parameter distribution A = A(r) in 
the limit r«T c . Although the structures of mutliquan- 
tum vortices have been calculated self-consistently in the 
framework of Bogolubov-de Gennes theory the vortex 
core anomalies have not been discussed ye t 20 ' 21 . More- 
over multiple anomalous branches of electronic spectrum 
have been shown to produce complicated patterns in 
the LDOS distributions investigated in the framework of 
Bogolubov- dc Gcnncs theor y 20 ! 21 . Here we employ an al- 
ternative approach of quasiclassical Eilcnbcrger theory^ 2 - 
to check the predictions of vortex core anomalies and the 
LDOS interference patterns in multiquantum vortices in 
clean superconductors. 

Notwithstanding the interesting physics taking place 
in the clean regime the experimental realization of STM 
measurements of multiquantum vortex states was im- 
plemented on Pb superconductor^*^ with short mean 
free path being much smaller than the coherence length. 
This dirty superconductor is more adequately described 
within the diffusive approximation of the electronic mo- 
tion resulting in the Usadel equations for the electronic 
propagators and the superconducting order parameter—. 
Singly quantized vortex states in dirty superconductors 
was investigated in detail^ and were shown to lack the 
low temperature singularity of the A(r) distribution be- 
ing smoothed out by the impurity scattering of quasi- 
particle states. Moreover the LDOS distribution inside 
vortex core does not feature zero bias anomaly since the 
spectral weight of bound electronic states is distributed 
smoothly between all energy scales up to the bulk energy 
gap Ao. On the other hand the multiquantum vortex 
states have not been investigated in the framework of the 
Usadel theory nor the LDOS distributions around mul- 
tiquantum vortices in dirty superconductors have been 
ever calculated. 

It is the goal of the present paper to study both the 
peculiarities of the multiquantum vortex structures espe- 
cially at low temperatures and the distinctive features of 
the electronic LDOS near the vortices which would allow 
unambiguous identification of giant vortices both in clean 
and dirty regimes. This paper is organized as follows. In 
Sec. [II] we give an overview of the theoretical framework 
namely the quasiclassical Eilenberger theory in clean su- 
perconductors and Usadel equation in the dirty regime. 
We discuss the results of self-consistent calculations of 
the order parameter distributions for multiquantum vor- 



tex configurations in Sec. IHII and address the LDOS pro- 
files in Sec. IIV1 We give our conclusions in Sec. [V] 



II. THEORETICAL FRAMEWORK 

A. Clean limit: Eilenberger formalism 

Within quasiclassical approximation the band param- 
eters characterizing the Fermi surface is the Fermi ve- 
locity Vf and the density of states No. We normalize 
the energies to the critical temperature T c and length to 
£c = ^Vf/T c . The system of Eilenberger equations for 
the quasiclassical propagators reads 



n p (V + iA)f + 2wf - 2Ag = 0, 
n p (V - iA) /+ - 2ujf+ + 2A*g = 0. 



(4) 



Here A is a vector potential of magnetic field, the vector 
rip parameterizes the Fermi surface and u is a real quan- 
tity which should be taken at the discrete points of Mat- 
subara frequencies u n = (2n + l)7rT determined by the 
temperature T. The quasiclassical propagators obey nor- 
malization condition g 2 + ff + = 1. The self-consistency 
equation for the gap is 



A(r) 



N d 

„=n JFS 



f(u n ,r,n p )d 2 S p . 



(5) 



where A is coupling constant, Sf is a Fermi surface area 
and the integration is performed over the Fermi sur- 
face. Hereafter to simplify the calculations we assume 
the Fermi surface to be cylindrical and parameterized 
by the angle 6 P so that n p = (cos8 p ,sm6 p ). In Eq.® 
Nd(T) = Wd/(2vrT) is a cutoff at the Debye energy ojd 
which is expressed through physical parameter T c and A 
as follows 



N d (T c ) 

E ■ 



A 



1/2 



= 1. 



(6) 



The LDOS is expressed through the analytical continu- 
ation of quasiclassical Green's function to the real fre- 
quencies 



Sp 1 I Re[g{uo = -iE + 0, r, rip)]^^. (7) 

JFS 



N{v) = N 



Assuming the vortex line to be oriented along the z axis 
we choose the following ansatz of the superconducting 
order parameter corresponding to axially symmetric vor- 
tex bearing M quanta of vorticity A(x,y) = |A|(r)e lMv 
where r = \/ x 2 + y 2 is the distance from the vortex cen- 
ter, ip = arctan(y/x) is the polar angle. Below we ne- 
glect the influence of the magnetic field on the vortex 
structure which is justified for superconductors with large 
Ginzburg-Landau parameter. 
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B. Dirty limit: Usadel equations 

In the presence of impurity scattering the Eilenberger 
Eqs. contain an additional diagonal self-energy 

term^ 2 -. When the scattering rate exceeds the corre- 
sponding energy gap (dirty limit) the Eilenberger theory 
allows for significant simplification. In this case the qua- 
siclassical Usadel equations 2 ^ are applicable. The struc- 
ture of singly-quantized vortices with M = ±1 in dirty 
superconductors was studied extensively in the frame- 
work of the Usadel equations^ 

ujF - [G(V - iA) 2 F - FV 2 G] = AG (8) 

where G and F are normal and anomalous quasiclassical 
Green's functions averaged over the Fermi surface sat- 
isfying the normalization condition G 2 + F*F = 1. To 
facilitate the analysis, we introduce reduced variables: we 
use T c as a unit of energy and £,d = \fT>/2T c where V 
is a diffusion constant as a unit of length. The Usadel 
equation is to be supplemented with the self-consistency 
equation for the order parameter 

N d 

A(r) =27rTA^ J FV„,r). (9) 

We again neglect the influence of the magnetic field on 
the vortex structure. It is convenient to introduce the 
vector potential in Eq. ([5} corresponding to a pure gauge 
field which removes the phase of the order parameter 

A = M^. (10) 

Using 8— parametrization^ (F = sin#, G = cos 9) the 
Usadel equations and self-consistency conditions can be 
rewritten in the form 

I A. (r4-o) -^sin(20) + (Acos0-o;siii0) = 0. 
r dr \ dr J 2r 2 

(11) 

Local density of states (LDOS) N(E,r), which is ac- 
cessible in tunneling experiments, can be obtained from 
9{uj, r) using analytic continuation 

N(E,r) = Rc[cos6»(w -> -iE + S,r)} (12) 

III. ORDER PARAMETER STRUCTURES OF 
MULTIQUANTUM VORTICES 

To determine the behavior of gap functions A = A(r) 
in multiquantum vortices we solved numerically the sets 
of Eilenberger Eqs. (|4I5[) and Usadel Eqs. (|9llip which 
describe the clean and dirty regimes correspondingly. At 
first let us consider the clean regime. The order pa- 
rameter profiles in vortices with winding numbers M = 
1, 2, 3, 4 are shown in Fig[TJa,b,c,d) for the temper- 
atures T/T c = 0.1; 0.5; 0.9. One can see that at ele- 
vated temperatures T = 0.9T C (red dashed curves) and 



T = 0.5T C (green dash-dotted curve) the order parame- 
ter follows Ginzburg-Landau asymptotic A(r) ~ r AI at 
small r. 

At low temperature T = 0.1T C the order parameter dis- 
tribution inside vortex core is drastically different from 
the Ginzburg-Landau behaviour as shown by blue solid 
lines in Fig.([T|). In particular the singly quantized vortex 
in FigQJa) features the Kramer-Pesch effect^ when the 
order parameter slope at r — grows as dA/dr ~ 1/T 
when T — > 0. In case of mutiquantum vortices with 
A I > 1 the anomalous massless branches of electronic 
spectrum (J3j> produce anomalies in the multiquantum 
vortex core structures 2 ^. To observe the vortex core 
anomalies we plot in Fig.@ the derivatives dA(r)/dr of 
the order parameter profiles obtained self consistently for 
the vortex winding numbers M = 1, 2, 3, 4. In accordance 
with the analytical consideration the vortex core anoma- 
lies result in the singular behavior of dA(r)/dr at low 
temperatures. We find that at T = 0.1T C in multiquan- 
tum vortices with M > 1 the calculated dependencies 
dA(r)/dr have sharp maxima at finite r ^ 0. Accord- 
ing to the analytical predictions these maxima originate 
from the square root singularity of the order parameter 
which is produced by the contribution of the anomalous 
energy branches of electronic spectrum 2 ^. 

In general for higher values of winding numbers M > 1 
in the limit T —> one should have M/2 singularities of 
dA{r)/dr for even M and (M+ 1)/2 singularities for odd 
M. For the particular examples of M = 2,4 there are 
one and two peaks of dA/dr at T — 0.1T C shown by blue 
solid line in Fig(5](b,d). Vortices with odd winding num- 
bers should feature Kramer-Pesch effect analogous to the 
singly quantized vortex originating from the anomalous 
energy branch crossing the Fermi level at (i = in the 
Eq. ([3]) . We found that the order parameter of M = 3 vor- 
tex has linear asymptotic A(r) ~ r at small r shown in 
the FigQTc). The slope of this linear dependence should 
grow at decreasing temperature which is an analog of 
Kramer-Pesch effect. The detailed analysis of vortex core 
anomalies in multiquantum vortices and in particular of 
the Kramer-Pesch effect for odd winding numbers M will 
be published elsewhere. 

Next consider the case of dirty superconductor and cal- 
culate the core structures of multiquantum vortices. The 
results of calculation are shown in Fig.(j3]) for the wind- 
ing numbers M = 1,2,3,4 and temperatures T/T c = 
0.1; 0.5; 0.9. As expected the vortices in dirty regime 
do not feature singularities in the order parameter dis- 
tribution in contrast to the clean case considered above. 

The comparison of vortex core structures in clean and 
dirty superconductors at T/T c = 0.1 is presented in Fig|3] 
for the winding numbers M = 1,2,3,4. To demon- 
strate the difference between clean and dirty cases we 
plot the dependencies A = A(r) in logarithmic scale in 
FigsQJb,d) correspondingly In the dirty case the order 
parameter has Ginzburg-Landau power law asymptotic 
A(r) = ar M which takes place at r — > even at very low 
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(a) M=1 (b) M=2 




1 2 3 1 2 3 




1 2 3 1 2 3 



FIG. 1: The distribution of the order parameter around vor- 
tex cores in clean superconductor at different temperatures. 
The panels (a,b,c,d) correspond to the winding numbers M — 
1, 2, 3, 4. Blue solid, green dash-dotted and red dashed lines 
correspond to the temperature T/T c = 0.1; 0.5; 0.9. 

temperatures TcT c . In Fig|4ja,b) the low-temperature 
behavior A(r) in the clean case is drastically different 
from Ginzdurg-Landau regime. In particular the multi- 
quantum vortex with M — 3 shown by blue dash-dotted 
line in Fig|4^ has linear asymptotic at r = 0. The slope 
of linear asymptotic for M = 3 should grow with decreas- 
ing temperature featuring an analog of Kramer-Pesch ef- 
fect for multiquantum vortices. Furthermore the order 
parameter in M = 4 vortex shown by red dashed line 
in FigH^i is almost zero at finite region r < R c where 
Rc ~ £c/2- This behavior is caused by the dominat- 
ing contribution of the electronic states corresponding 
to anomalous branches to the superconducting or- 
der parameter at r < R c . Thus contribution is zero at 
r < min((ioi, [iQzj/kF in the limit T —> 021. Thus the 
multiquntum vortices with even winding numbers M are 
well described by the step-wise vortex core model used 
previously for the analytical analysis of the vortex core 
spectrum^. 



IV. LDOS PROFILES OF MULTIQUANTUM 
VORTICES 

Having in hand the order parameter structures calcu- 
lated self-consistently is Sec. (|III[) we calculate the LDOS 
distributions formed by the electronic states localized at 
the vortex core. We start with the case of clean super- 
conductor which is known to demonstrate peculiar pro- 



fa) M=1 (b) M=2 




1 2 3 1 2 3 



FIG. 2: The vortex core anomaly revealed at the sharp 
peak of radial dependence of the order parameter profile 
derivative dA(r)/dr normalized to the value ^c/T c around 
vortex cores in clean superconductor at different tempera- 
tures. The panels (a,b,c,d) correspond to the winding num- 
bers M = 1,2,3,4 correspondingly. Blue solid, green dash- 
dotted and red dashed lines correspond to the temperature 
T/T c = 0.1; 0.5; 0.9. 



files of LDOS originating from multiple anomalous en- 
ergy branches of localized electron o 20 ' 21 . Here we cal- 
culate the LDOS distributions for the winding numbers 
M = 1, 2, 3, 4 shown in Fig.©. The LDOS plots are sim- 
ilar to that obtained in the framework of Bogolubov- de 
Gennes theor y 20 ' 21 . 

Introducing a polar coordinate system (r, ip) and defin- 
ing the z projection of quasiparticlc angular momen- 
tum through the impact parameter of quasiclassical 
trajectory 2 ^ /x = — [r, kp] • zo the LDOS inside the 
singly quantized vortex core can be found with the help 
of Eq.@ as follows: N{E,r) = (k F /2n£ c ) 5 [ E ~ 
ojfcprsin(<^ — 8 p )]d6 p . Here we evaluate the LDOS sum- 
ming up over the quasiparticle states at the trajectories 
characterized by the direction of the quasiparticle lin- 
ear momentum kp = kp(cos 9 P , sin 9 P ). This expression 
yields a singularity of zero energy LDOS^Ir— : N(E, r) — 
l/(2^ v ^ 2 ^7 2 e c ) « iVoCc/v/^7 2 , where N = 
(l/2ir)m/H 2 is a normal metal LDOS and ro = E/iukp). 
Thus the LDOS profile of singly quantized vortex has the 
ring form with the radius r being a function of energy. 
The dependence N = N(E,r) is shown in Fig. ©a for a 
singly quantized vortex. 

In multiquantum vortices the spectrum of low energy 
states contains several anomalous branches which 
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(a) M=1 



(b) M=2 



(a) Clean limit 



(b) Clean limit 




FIG. 3: The distribution of the order parameter around 
vortex core in dirty superconductor at different tempera- 
tures. The panels (a,b,c,d) correspond to the winding num- 
bers M = 1,2,3,4 correspondingly. Blue solid, green dash- 
dotted and red dashed lines correspond to the temperature 
T/T c = 0.1; 0.5; 0.9. 



intersect the Fermi level and contribute to the low- 
energy DOS. The LDOS profile corresponding to the 
spectrum ([3]) consists of a set of axially symmetric ring 
structures^—. Note that for an even winding number 
the anomalous branch crossing the Fermi level at /i = 
(i.e. at zero impact parameter) is absent and, as a re- 
sult, the LDOS peak at the vortex center disappears. 
Using the same procedure as for the singly quantized 
vortices and the spectrum ^ we obtain the LDOS in 

" 2 i£c) where 



the form N(E,r) = £ -=i V( 27I "w 4 ypr" - 
Tai = [A'oi + E/uji\/kp. At E = the spectrum is sym- 
metric so that the LDOS profile has M/2 peaks for even 
M and (M + l)/2 for odd M. At E ± the degeneracy 
is removed and each peak splits by two as can be seen 
from the LDOS plots in Fig.©. 

Smearing of energy levels due to scattering effects leads 
to a reduction of LDOS peak at the vortex center. How- 
ever, the LDOS peak survives even in "dirty" limit when 
a mean free path is smaller than a coherence length I < £. 
To find the form of LDOS peak at the vortex core we 
consider the dirty case described by Usadel Eq. (fTTj) . The 
LDOS distributions around multiquantum vortices cal- 
culated according to Eqs. (|llll2[) are shown in Fig.®. 
The profiles N = N(r) in multiquantum vortices M > 1 
feature plateau at zero energy level. This is in high con- 
tract to the case of singly quantized vortex M = 1. In 
fact these plateau occur due to the very slow spatial dc- 




FIG. 4: The distribution of the order parameter around mul- 
tiquantum vortex core at T/T c = 0.1 in (a) dirty supercon- 
ductor and (b) clean superconductor. Black solid, green dot- 
ted, blue dash-dotted and red dashed lines correspond to the 
winding numbers M = 1, 2, 3, 4. 



(a) M=1 




(b) M=2 




(d) M=4 




FIG. 5: The distribution of the LDOS around vortex cores 
at T/T c = 0.1 in clean superconductor as function of energy 
and distance from the vortex core iV = N(r,E). The panels 
(a,b,c,d) correspond to the values of vorticity M = 1,2,3,4 
correspondingly. 



pendence of SN(r) = 1 — N(r)/N at small r which can 
deduced directly from Eqs. (|llll2l) . Indeed linearizing the 
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Eg .([IT]) for 
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we obtain 


'1 


(I 


( d 9] AP 
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dr 


\ dr J r 2 



A(r) 



= 



(13) 



which define the asymptotic 9(r) = ar . Next the Eq. 
(fT2]) yields the LDOS SN = 9 2 /2 = a 2 r 2M /2. This an- 
alytical asymptotic perfectly agrees with the numerical 
results which can be seen from the logarithmic scale plot 
of N(r) in Fig.© b. An interesting feature of such LDOS 
plateau is that they survive at the distances compared to 
the size of the multiquantum vortex core which is much 
larger than the coherence length. That is we find that 
the size of the plateau shown in Fig.© is approximately 
with good accuracy R p = M£p/2 for M > 1. 



(a) M=1 



(b) M=2 




0.5 

E/A 15 10 5 
(C) M=3 





FIG. 6: The distribution of the LDOS around vortex cores 
at T/T c = 0.1 in dirty superconductor as function of energy 
and distance from the vortex core N = N(r,E). The panels 
(a,b,c,d) correspond to the values of winding number M — 
1, 2, 3, 4 correspondingly. 



V. CONCLUSION 

To summarize we have calculated self-consistently in 
the framework of quasiclassical Eilcnberger theory the 
order parameter structures of multiquantum vortices to- 
gether with the local density of electronic states both in 
clean and dirty superconductors. We have fond that at 
the temperatures near T c the order parameter profiles 
of vortices are qualitatively similar in clean and dirty 
regimes (compare the plots by red curves for T = 0.9T C 
in FigsfT] and [3]). In this temperature regime the or- 
der parameter asymptotic is determined by the power 




r/^ D 10 



^0 10 2 



FIG. 7: (a) The LDOS profile for zero energy E = around 
vortices at T/T c = 0.1 in dirty superconductor as function 
of the distance from the vortex center N = N(r). (b) The 
logarithmic plot of SN(r) = 1 — N(r)/No demonstrating the 
power law asymptotic SN(r) ~ r 2M at r — > 0. Black solid, 
green dotted, blue dash-dotted and red dashed lines corre- 
spond to the winding numbers M = 1, 2, 3, 4. 



law A(r) = ar M which is consistent with the result of 
Gizburg-Landau theory valid at \T/T C — 1| <C 1. 

On the other hand in the low temperature limit 
T = 0.1T C vortices in clean superconductor demon- 
strate the anomalies in the order parameter distribu- 
tion - the singularities of the derivative dA/dr predicted 
in Ref.(d3). Such singularities occur due to the contri- 
bution of anomalous electronic spectral branches to the 
superconducting order parameter. The singular behav- 
ior of dA/dr in multiquantum vortices is analogous to 
the Kramer-Pesch effect^ taking place for singly quan- 
tized vortex M = 1 which has steep order parameter 
slope dA/dr(r = 0) ~ l/T at T -> (compare the blue 
solid curves which correspond to T/T c = 0.1 in FigsQJ, 
and [SJd). In dirty superconductors the order parame- 
ter asymptotic at the vortex core obeys the Ginzburg- 
Landau power law behavior even at low temperature 
T = 0.1 T c which is clearly demonstrated in logarithmic 
scale plots in FigHJi. 

In the framework of quasiclassical theory we calculated 
the LDOS distributions in muliquantum vortices with 
winding numbers M = 1,2,3,4. The LDOS profiles in 
the clean regime are similar to that obtained previously 
with the help of Bogolubov-de Gennes theor y 20 i 21 . Most 
importantly we determined the LDOS profiles in dirty 
regime which directly correspond to the modern exper- 
iments on scanning tunneling microscopy of multiquan- 
tum vortices in mesoscopic superconductors. The asymp- 
totic of LDOS profile is shown to be N(r)/N = l — ar 2M 
which holds with good accuracy at r < M^d/2. Thus for 
the values of M > 2 the LDOS profile is almost flat at 
the sizable region near the vortex center r < M£d/2 (see 
Fig II]). Such LDOS plateau can be considered as a hall- 
mark of multiquantum vortices in STM experiments in 
dirty superconductor o 9 ' 10 . 
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